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Based on first-principles calculations, we predict a substantial increase in the optical di¬ 
electric function of LiAsSe2 under pressure. We find that the optical dielectric constant 
is enhanced threefold under volume compression. This enhancement is mainly due to the 
dimerization strength reduction of the one-dimensional (ID) As-Se chains in LiAsSe2, which 
significantly alters the wavefunction phase mismatch between two neighboring chains and 
changes the transition intensity. By developing a tight-binding model of the interacting ID 
chains, the essential features of the low-energy electronic structure of LiAsSe2 are captured. 

Our findings are important for understanding the fundamental physics of LiAsSe2 and pro¬ 
vide a feasible way to enhance the material optical response that can be applied to light 
harvesting for energy applications. 
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I. INTRODUCTION 


The dielectric response, as a fundamental physical property of materials, describes how ma¬ 
terials respond to an external electric held. In semiconductors, when the applied electric held 
frequency is in the range of visible light, the photon excitation of electronic inter-band transitions 
dominates the total dielectric response, which is described by the optical dielectric function. The 
optical dielectric function is strongly related to other optical properties of the material, including 
light absorption, refraction and non-linear optical responses. Therefore, the enhancement and tun- 
ability of the optical dielectric function of a material are signihcantly important in various areas, 
such as solar cell, optical devices and sensors. A great deal of research has been done to increase 
the material optical dielectric response. In particular, defects, material doping and surface plas- 
mon induced by metallic nanoparticles have been widely used to increase the optical absorption 
in semiconductorsi”— . Whereas most of the previous methods rely on the assistance of another 
material, the intrinsic bulk dielectric response enhancement of the light absorber is less studied. 

Alkali-metal chalcogenides such as KPSeg, K2P2Se6, LiAsSe2, LiAsS2 and NaAsSe2 have been 
synthesized, and their band gaps lie in the visible light region^. Since they have spontaneous polar¬ 
ization, these materials are potential candidates to show the bulk photovoltaic effect^. Moreover, 
strong optical second-harmonic generation susceptibility has been observed experimentally and the¬ 
oretically^*^. However, the effect of structural distortion on their linear optical responses has not 
been studied, and the structure-property-optical performance relationship is still unclear. In this 
paper, by using a first-principles method, we show that the optical dielectric constant of LiAsSe2 
increases threefold by volume compression. More interestingly. As and Se atoms in LiAsSe2 form 
weakly interacting quasi-one-dimensional atomic chains, of which the dimerization strength can be 
tuned by volume compression. Atomic chains have attracted a great deal of interest, due to their 
one-dimensional nature giving rise to exotic phenomena such as conductivity^*^, metal-insulator 
transitionii, and topological phase o^^*^^ . Herein, their important roles in light absorption are em¬ 
phasized. As illustrated by a tight-binding model, the dimerization strength is strongly coupled to 
the relative phases of the gap state wavefunctions between the two neighboring chains. By reducing 
the wavefunction phase mismatch between the chains, the magnitude of transition intensity for the 
transitions near the band edges increase significantly, giving rise to substantial optical dielectric 
function enhancement. 


3 


II. COMPUTATIONAL DETAILS 


Figuresand b show the experimental structure (ES) of LiAsSe2^. The polar phase of LiAsSe2 
has the Cc space group with the glide plane perpendicular to the lattice vector h. The polarization 
induced by ionic displacement lies in the a-c planed. As shown, the As and Se atoms form distorted 
quasi-one-dimensional atomic chains along the h direction. This chain and its neighboring chains 
form a two-dimensional chain plane (illustrated as the grey plane), and these chain planes are 
separated by Li-Se planes (light purple plane). In the ES, this As-Se chain dimerizes, creating 
alternating As-Se bonds with two different bond lengths. 
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FIG. 1 . a) The unit cell of LiAsSe2. The lines between As and Se atoms indicate the quasi-one-dimensional 
chains. The chain with its neighbor chains form a chain plane (grey color plane). These parallel chain 
planes are separated by the Li-Se plane (light purple plane) in the middle, b) Side view of the experimental 
structure (ES). c) Side view of the compressed structure (CS). The differences between the ES and the CS 
are mainly ion motions in the b direction. As illustrated by the bond lengths between two neighboring As-Se 
bonds, ES shows stronger dimerization strength along the chain than CS. 


The plane-wave density functional theory (DFT) package QUANTUM-ESPRESSO was used 
to perform structural relaxations and electronic structure calculations, with the Perdew-Burke- 
Ernzerhof (PBE) generalized gradient approximation exchange-correlation functional^^. Norm- 
conserving, designed non-local pseudopotentials were generated with the OPIUM package^*^. A 
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plane-wave cutoff energy of 50 Ry was sufficient to converge the total energy with the fc-point 
sampling on a 4x8x8 grid. The structure relaxed with the PBE functional underestimates the 
dimerization along the chain, and it does not match with the ES. By using the GGA -|- U method 
with effective Hubbard Ueff = 7.5 eV on the As 4 p orbitals, the relaxed structure matches the ES 
very well. Adding U on p orbitals to get the correct structure is not rare, as the large self-interaction 
error originating from s or p orbitals may partially be corrected by the DFT-|-?7 metho d^^'^^ . The 
DFT calculated band gap is 0.8 eV, which underestimates the experimentally measured 1.1 eV— . 
With the converged charge density, the wavefunctions used for the dielectric function calculations 
are obtained from non-self-consistent calculations performed on a denser fc-point grid with a suffi¬ 
cient number of empty bands. By using the long wavelength approximation and the single particle 
approximation, the imaginary part of the optical dielectric function is calculated as Eq.([T|), 
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where oj is the light frequency; i is the Gartesian coordinate; k is the Bloch wave vector; c, v denote 
the conduction and valence band with energy hucfv The real part of the dielectric function, ei, 
can be calculated from the Kramers-Kronig relation. 


III. RESULTS AND DISCUSSION 

The compressed structure (GS), with much weaker dimerization strength of the atomic chains 
(Fig. [It), is obtained by compressing all the lattice vectors by 3%, followed by the relaxation of 
the internal atomic positions. This volume compression of LiAsSe 2 strongly enhances its optical 
dielectric response as shown in the calculated optical dielectric functions of the ES and GS (Fig. [2t). 
Figures |2]3, c and d illustrate the calculated joint density of states (JDOS), refractive index and 
absorption spectrum along the b direction as a function of the photon energy for the ES and GS, 
respectively. As shown from the spectrum, the GS shows much higher linear optical responses 
near the band gap than the ES. In particular, the optical dielectric constant of the GS increases 
to more than three times its original value (Fig. [2^). The imaginary part of the optical dielectric 
function, describing the real electronic inter-band transitions, also shows great enhancement under 
compression. As expressed in Equation ([T]), the imaginary part of the dielectric function is the 
product of JDOS He i; k '^(^c,k “ ^v,k — and the transition intensity \{c,k.\pi |u,k)p. However, 
we find that the JDOS contribution to the enhancement is negligible. As shown in Fig. [2]3, in 



5 



FIG. 2. a) LiAsSe 2 optical dielectric (e) function spectrum of the ES and CS as a function of photon energy 
along b. ei is the real part of the dielectric response spectrum, and £2 is the imaginary part, b) Joint density 
of states for the two structures. Owing to the different band gaps of the ES and CS (0.2 eV difference), 
the inset graph shows the shifted-CS (shifting the spectrum by 0.2 eV) and ES JDOS spectrums in order to 
compare with the same band gaps, c) Refractive index (n) spectrum along b. d) Absorption coefficient (a) 
spectrum b. 


the energy range 0 < huj < 2 eV where the imaginary part €2 shows substantial enhancement, 
the calculated JDOS for the ES and CS have very similar magnitude. Therefore, this dielectric 
function enhancement mainly comes from the increase of transition intensity by compression. 

To further resolve the origin of the dielectric enhancement by pressure, we present the distribu¬ 
tion of the transition intensity as a function of k in momentum space. Figure [3] shows the transition 
intensity distributions in the Brillouin zone (BZ), with the transitions between the valence band 
maximum (VBM) to the conduction band mimimum (CBM) within the the energy range of 0-2 
eV, since these transitions are dominant in the dielectric function enhancement. As displayed in 
Figure El the /c-resolved distributions show distinct patterns in addition to their overall differ¬ 
ences in the corresponding dielectric constants. For the ES, most of the k points have similar yet 
low magnitude of transition intensities. However, for the CS, the fc-resolved transition intensity 
shows significant changes, with the high magnitude k points mostly distributed on a thin plane 
perpendicular to the reciprocal lattice vector kb. The k points contributing the highest transition 
intensities are broadly located in the middle region in this plane. Along the ka and kc directions 
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of this plane, the transition intensity changes slowly with respect to wavevectors, indicating the 
weak bonding character. This can be attributed to the weak As/Se-Se/Li inter-planar and As-Se 
inter-chain interactions. However, the magnitude of transition intensity shows rapid change along 
the kb direction, as illustrated by the transition intensity profile along this direction (Fig. [ 6 ) 3 ). 
This strong fc-dependent transition intensity distribution reveals the strong covalent bonds char¬ 
acter along the chain direction. The highly inhomogeneous distribution of the transition intensity 
can be considered as an indication of the quasi one-dimensional nature of the system near the 
low-energy spectrum, steming from the dimerization changes of the As and Se atoms. Therefore, 
investigating the electronic structure of the chains is essential to further understand the origin of 
the dielectric response enhancement. 



FIG. 3. Distribution of |(V’w,k| P \'4’c,'k)f jV (eV/A^) in the Brillouin zone (BZ) extracted from DFT calcu¬ 
lation of LiAsSe 2 for a) the ES and b) CS. V is the volume of the unit cell. For simplicity, the primitive 
BZ is illustrated as an orthogonal box with reciprocal lattice vectors ka, kb, and kc along the three edges of 
the box. The transitions with transition energy less than 2 eV are plotted, as this energy region shows the 
greatest dielectric function enhancement. The detailed transition intensity profiles along the black lines in 
the hgures for the ES and CS are shown in Fig.[B}3. 


Figure S] shows the DFT band structures plotted along F-Y. Under compression, most bands 
along the kb direction show relatively small changes, except the bands near the band edges. The CS 
shows strong dispersion near its optical gap at (0.0,0.42,0.0) (fractional coordinate), as its CBM 
shows a “dip” while the VBM shows a “bump”. The ES has its band gap shifted towards the BZ 
boundary at (0.0,0.46,0.0). Comparing to those of the CS, the bands of the ES near the band gap 
shows much less dispersion, and the dip and bump features become less obvious. Besides the band 
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dispersion change, the band gap shows noticeable change from 0.80 eV (ES) to 0.62 eV (CS). More 
importantly, we find that the inter-band transition between the band edges in the CS provides the 
highest transition intensity magnitude, but this corresponding value in the ES is very low. In order 
to understand the bonding characters of these states which give the highest transition intensity, the 
charge density isosurfaces of the VBM and CBM are plotted in FigureSl Unexpectedly, both the ES 
and CS show quite similar charge density distributions, with non-bonding Se p orbital character as 
VBMs and non-bonding As p as CBMs, suggesting that the atomic orbital overlaps cannot explain 
such large dielectric enhancement by compression due to their similar charge densities. Rather, 
we find that the dimerization change induced by the compression can strongly alter the phase of 
wavefunctions so as to vary transition intensity magnitude significantly, as we will discuss below. 
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FIG. 4. The band structure of LiAsSe 2 from F to Y (0, n/b, 0) along kb, and the charge density isosurfaces 
of the conduction band minimum (CBM) and valence band maximum (VBM) states indicated by the blue 
squares in the band structures for a) ES and b) CS. 

To demonstrate the significant influence of wavefunction phase change on the optical response 
enhancement, we construct a two-dimensional (2D) tight-binding (TB) model with interacting 
atomic chains illustrated in Figure [5^. The TB model comprises four orbitals (j =1, 2, 3, 4) in 
a square lattice with lattice constant a and periodic boundary conditions along b and c to model 
a chain plane in LiAsSe 2 . Owing to the weak interaction between the chain plane and the Se-Li 
plane, the inter-planar interaction along the a direction is not considered. As shown in the charge 
density distributions (Fig. [3|), the p orbitals from the As and Se atoms form u-type covalent bonds 
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FIG. 5. a) 2D TB model for weakly interacting As-Se chains (the inset graph shows the chains in LiAsSe 2 ). 
The dashed lines indicate the chain-chain interaction connecting the As-Se chains (solid lines). tzL6t denotes 
the hopping strength, b) The Brillouin zone of the 2D model. The band structure (graph c) is plotted along 
the thick blue line, c) The band structure calculated from the 2D TB model along the chain propagation 
direction under different dimerization strengths {Sti/t 2 with t 2 fixed). 

along the chain (along the b direction). The dimerized hopping strength is denoted as ti ± 5ti to 
describe the alternating As-Se bond lengths. Across the chains (c-direction), vr-bonding between 
the p orbitals forms, where the corresponding hopping interaction is denoted as t 2 ± dt 2 - We 
find that this inter-chain interaction is of crucial importance in reproducing the correct DFT band 
structure, although these interactions are weak relative to the intra-chain interaction, thus assuming 
1 ^ 2 ! < |G|- The onsite energies of As and Se orbitals are set to Eq + 6E and Eq — SE, respectively. 
The onsite energies and hopping strengths of the TB Hamiltonian are tuned to reproduce the DFT 
band structure. 

By solving the TB model numerically, we obtain the band structures in Fig. [dh plotted along 
the chain propagation direction indicated by the blue line in Fig. [5)3. We also calculate the band 
structures by gradually reducing the dimerization strength (decreasing 6ti) with t 2 fixed, and find 
that the band gap position shifts away from the BZ boundary. More interestingly, the dispersion 
of the band edges are significantly enhanced when decreasing 5ti. This feature becomes clearer by 
calculating the fc-resolved transition intensity using the TB wavefunctions. 

As derived in the Appendix, the transition intensity (X) is expressed as I (k) = | (k) H (k)|^, 

with ji C'J*kC'jk’ summing over the contributions of the wavefunction coefficients to the 

transition intensity. In this case, H accounts for the contribution generated when constructing 
the TB basis set from the localized atomic orbitals, with its relative value only determined by 
the wavevector without solving the Hamiltonian. Shown in Fig. [6] is the calculated I along the 
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same A:-path as used in the TB band structure (Fig. [ 5 ) 3 ). By reducing the dimerization strength of 
the atomic chains (reducing 5 ti/t 2 ), the transition intensities for the band edge states and nearby 
increase significantly, which agrees well with the DFT transition intensity trend under compression 
shown in Fig. [ 6 ) 3 . Additionally, by plotting W which is contributed only from the wavefunction as 
we are interested, it is clear that it shows exactly the same trend as I, demonstrating the significant 
role of wavefunctions in the enhancement of the transition intensity under pressure. 

The low-energy k ■ p effective theory provides simpler and more explicit band structure and 
wavefunction expression. The Hamiltonian H (k) is further expanded in the vicinity of the BZ 
boundary as ^^(k) = H(K) -|- (k — K)ff'(K) with k = K -|- (g, 0), K = ( 7 r/ 6 ,0). From the k ■ p 
Hamiltonian, the energies for the valence band (FI-) and the conduction band (Fi_|_) near the BZ 
boundary are obtained as: 


E± (q) = ±\j5E‘^-2^JulVL (q) + AStj + H (q) (2) 

H (q) = + {qat{f (3) 

When \Sti\ > \t 2 \, the band gap is at the BZ boundary (g = 0). By decreasing the dimerization 
strength such as \5ti\ < \t 2 \, the band gap wavevector {q{Eg)) is 2 ^/t 2 — Stf/ (ati). This change of 
band gap position as a function of the dimerization strength { 6 ti/t 2 ) agrees with our DFT band 
structures of LiAsSe 2 . In the ES, the strong dimerization between the As and Se atoms moves the 
band gap close to the BZ boundary. In the CS, the reduced dimerization due to the compression 
shifts the band gap away from the boundary, giving rise to the strong dispersion for the states near 
the gap. 

Within this low energy theory, the phase relationships of the wavefunctions are further explored 
by evaluating the analytical expression of wavefunctions for the band edge states. When the band 
gap is not at the BZ boundary (|(5ti| < |t 2 |) j the wavefunctions of the gap states have simple 
forms; 


V’VBM = l/v^ (^0, e*^, 1, 0^ , (4) 

V-CBM = l/\/2 ( 1 , 0 , 0 , e'*^) (5) 


where the wavefunctions are written with the TB basis of the four orbitals: XAs,j=i5 XSe,j=2, XSe,j=3, 
and XAs,j =4 (Fig. [5^). Due to the simple form of the wavefunctions, we use these two states to show 
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FIG. 6. a) Calculated transition intensities I(k) and (k)|^ (k) = Yj,,' C'J'tkC'j.k) from 2D TB 

model. The x axis is the wavevector along the chain direction, b) Transition intensities (eV/A^) extracted 
from DFT calculation of transition intensity for the ES and CS. They are plotted along the chain direction 
as indicated by the short black lines in Fig. [31 


the effect of phases of the chains. Here, 9 = arcsin (|(5ti/t2|) indicates the dimerization strength 
of the atomic chains. From the wavefunction expression, it is clear that the VBM and CBM are 
always non-bonding states without mixing of the As and Se orbitals, which is also observed in the 
DFT calculation. More interestingly, 0 controls the phase mismatch between the wavefunctions 
of the two chains in the chain plane. For example, for the CBM wavefunction, when 0 = 0, 
the orbitals on XAs,j=i and XAs,j=A are populated in the same phase, while, with nonzero 0/0, 
XAs,j=i on one chain and XAs,j=4 on the neighboring chain have the phase difference of e*® between 
the corresponding wavefunction coefficients. Hence, the application of the hydrostatic stress to 
LiAsSe 2 reduces 0, enabling wavefunction phase matching between the two neighboring atomic 
chains, which is essential to the enhancement of dielectric responses. 


Using this simple form of wavefunctions, the band edge transition intensity is evaluated as 
X (X |e*® + e“®®|^ = cos^ (0) = 1 — The first exponential term e*® originates from one 

chain and the other term from the neighboring chain. In this material, the finite dimerization 
strength of the two neighboring chains have opposite effects to their contributions to the transition 
intensity. Therefore, without varying the overlaps between the atomic orbitals, the phase change 
of the wavefunction induced by structural change alters the overall dielectric function significantly. 
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IV. CONCLUSION 

In summary, by using a first-principles method, we have shown that volume compression can 
significantly enhance the optical dielectric function and the dielectric constant by factor of three 
in LiAsSe 2 . This material is essentially a network of As-Se ID atomic chains with the dimer¬ 
ization strength tunable by compression. The enhancement of the transition intensity near band 
edges is the main reason of the overall dielectric function improvement. A 2D tight-binding model 
with weakly interacting atomic chains is developed to explore the relation of dimerization strength 
and transition intensity. When the dimerization is strong, the wavefunctions of the two neighbor¬ 
ing chains have significant phase mismatch, providing destructive interference that reduces to the 
dielectric function. By reducing this wavefunction phase mismatch via compression, the collec¬ 
tive contributions from the chains dramatically enhance the overall dielectric response and light 
absorption. Our results indicate that this material is suitable as the light absorber in the solar 
cell application. Furthermore, since the transition intensity is related to other optical processes 
such as second-harmonic generation and the non-linear optical effects, we expect that the volume 
compression can enhance their responses. 
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APPENDIX 

The Bloch wavefunction based on the TB orbitals is: 


V’n.k = (A.6) 

3 

Xj (j=l, 2, 3, and 4) is expanded as and (/)rj' is the localized atomic orbital 

centering at the position of R -|- Sj. 
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With the Bloch wavefunctions, the transition intensity is expressed as; 


X(k) = |(V'^,k|p|'0c,k)P 

2 

= EC7':kCy.kn,y',x(k) (A.7) 

jj' 

where (k) = (0_r j/| p \ 4 >o,j) with summation over nearest hopping neigh¬ 

bor unit cells denoted by R, which is only related to the wavevector, orbital position and the 
momentum matrix element between two localized atomic orbitals. 

The low energy Hamiltonian is written as: 


E 

i25ti — qati 

2t2 

0 

i25ti — qati 

-E 

0 

2t2 

2t2 

0 

-E 

i26ti — qati 

0 

2t2 

—i25ti — qati 

E 


with respect to the four orbitals shown in Figure [5l Based on this Hamiltonian, the band edge 
states can be solved as Equations 4 and 5. 

When calculating the transition intensity for band edge transitions, the transition intensity can 
be further simplified as: 


2 




|n(g)r 


= \W-’^q)\^\Uiq)\^ (A. 8 ) 

In this model, nj=o,j'=i = The transition intensity is only related to the wavefunction 

coefficient C and the wavefunction phase mismatch between two neighboring chains. 
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